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Abstract. This paper presents and tests a previously unrecognised mechanism for driving a repli- 
cating molecular system on the prebiotic earth. It is proposed that cell-free RNA replication in 
the primordial soup may have been driven by self-sustained oscillatory thermochemical reactions. 
To test this hypothesis a well-characterised hydrogen peroxide oscillator was chosen as the driver 
and complementary RNA strands with known association and melting kinetics were used as the 
substrate. An open flow system model for the self-consistent, coupled evolution of the temper- 
ature and concentrations in a simple autocatalytic scheme is solved numerically, and it is shown 
that thermochemical cycling drives replication of the RNA strands. For the (justifiably realistic) 
values of parameters chosen for the simulated example system, the mean amount of replicant pro- 
duced at steady state is 6.56 times the input amount, given a constant supply of substrate species. 
The spontaneous onset of sustained thermochemical oscillations via slowly drifting parameters is 
demonstrated, and a scheme is given for prebiotic production of complementary RNA strands on 
rock surfaces. 
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1. Introduction 

Molecular self -replication is a fundamental process of all living organisms, and its initiation 
in the environment provided by the early Earth must have been essential for the emergence of life 
itself. In modern cells DNA replication is accomplished isothermally at moderate temperatures 
by a working army of enzymes using the free energy of ATP hydrolysis, and the strongly driven 
nonequilibrium condition is maintained by respiration and ion transport across cell membranes. 
However in the primordial soup there were no cells and no protein enzymes. The 'RNA world' 
hypothesis holds that cell-free RNA communities grew on solid surfaces and replicated, before 
the evolution of DNA JT). What energy source may have driven RNA replication in such an 
environment? 

It does seem that thermal cycling would be required to drive RNA replication in the absence of 
cellular (or other) machinery, because heat is required to dissociate double-stranded or multiplex 
RNA and a cooler phase is necessary for replication and annealing, given a supply of substrate 
template. This fact is often overlooked in hypotheses about the origin of life, as pointed out by 
Kovac et al. [2J. Exponential replication of RNA duplexes, fed by hairpins derived from an alanine 
tRNA and driven by externally imposed thermal cycling between 10 and 40°C, was achieved by 
Krammer et al. [3J. They proposed that, in the primordial soup, thermal cycling may have been 
provided by laminar convection between hot and cold regions in millimetre-sized rock pores flU, 
and that the substrate oligonucleotides could be concentrated by thermophoretic trapping. 

It has been suggested by Matatov B51 that chemical energy also may have played a role in 
prebiotic evolution, because, in the laboratory, heat liberated in the decomposition of aqueous 
hydrogen peroxide (H2O2) supported formation of amino acids from simpler precursors. In this 
work we propose that a natural mechanism for driving self-sustained thermal cycling is a thermo- 
chemical oscillator, and demonstrate in a specific case that the temperature cycling for driving the 
amplification of template RNA in environments where early life may have evolved can be provided 
by a thermochemical oscillator driven by exothermic reactions of H 2 0 2 . 

Exothermic reactions of aqueous H2O2 are well-known to give rise to robust, self-sustained 
thermochemical oscillations with frequencies around 0.02-0.005 s^ 1 [|6]-|9l- As our specific ex- 
ample we have chosen, with some justification expressed below, the oxidation of thiosulfate ion 
(S 2 03~) by H 2 0 2 , which is fast and highly exothermic: 

Na 2 S20 3 + 2H 2 0 2 -> ^Na 2 S 3 0 6 + ^Na 2 S0 4 + 2H 2 0. (Rl) 

The thermokinetics and thermochemistry of (Rl) have been well-studied, and experimental data 
from the literature are collected in table Q] The reaction is first order in both reactants. 

Oscillatory thermoconversion is typical of highly energetic, thermally sensitive, liquids such 
many of the peroxides ifTOl . The physical basis is as follows: In such liquids (and also in some 
vapor and gas phase substances), which have high specific heat capacity, the heat of reaction can 
be absorbed by the many intermolecular vibrational modes and intramolecular rotational modes. 
So the temperature rises but slowly as reaction proceeds, until these modes become saturated 
and consequently the temperature spikes — the heat cannot be absorbed any more by quantized 
excitations. But locally the reactant is depleted, so the temperature falls to a minimum before 



2 



Downloaded from http://biorxiv.org/ on September 18, 2014 



Table 1. Experimental data for reaction (Rl), presented in chronological order. 



AQmol-h- 1 ) E (kJ/mol) -AH (kJ/mol) RefT 



6.85 x 10 11 76.52 573.0 d 

2.13 x 10 10 68.20 585.8 d 

7.33 xlO 11 78.24 594.1 (H 

1.63 xlO 10 68.12 612.5 Q 

6.8 x 10 11 77.0 not measured El 

2.00 x 10 10 68.20 586.2 

2.0 x 10 10 68.3 562.8 JH1 



reactant accumulation allows reaction to occur and heat to be released, and the cycle begins again. 
With a steady supply of reactant these self-sustained thermal oscillations can continue indefinitely. 

1.1. H2O2 and thiosulfate on the early earth 

One school of thought holds that the geochemical environment for the emergence of life was 
provided by submarine hydrothermal systems and hot springs [fTTl . Experiments reported by Fous- 
toukos et al. llT8l strongly support the conjecture that H 2 0 2 is produced near hydrothermal vents 
when oxygenated seawater mixes with vent fluid. 

Another source of H2O2 production in such an environment involves a surface reaction of 
pyrite (FeS2) with H 2 0 and in experiments 0.391-0.567 mM/m 2 was produced lfl9l l20l . The 
specific surface area of pyrite is 2.0-4.0 m 2 /g fl2TTl . so there is a very real possibility that high 
concentrations of H 2 0 2 could build up locally and be supplied at a constant rate for long periods 
of time. A very credible body of work (see references in [fT9l ) holds that the most primitive 
photo synthetic cells used H 2 0 2 as an electron donor, so it is reasonable to assume that H 2 0 2 was 
produced long before those organisms evolved. Moreover it is believed that pyrite was present on 
the early earth. On the modern earth, pyrite deposits associated with hydrothermal activity can 
reach thicknesses of tens to hundreds of meters and spread over thousands of square kilometers 
within the crust Il22l . 

Prebiotic production of H 2 0 2 also may have occurred by photochemical disproportionation of 
the superoxide radical (0 2 ) in sunlit waters to H 2 0 2 [|23l . 

Thiosulfate ion is observed to occur in hydrothermal waters of Yellowstone National Park ll24l . 
so it is reasonable to conjecture that it was present in hydrothermal environments on the early earth. 

Porous rocks around hydrothermal vents therefore could provide microenvironments where 
naturally self-sustaining thermochemical oscillations may be set up — namely, a localized, nonequi- 
librium forced flow system, and a supply of H 2 0 2 and thiosulfate ion. If a supply of RNA oligonu- 
cleotides, perhaps synthesized on and dislodged from pore surfaces 11251 . is added to this recipe 
for primordial soup, replication may be driven by thermochemical temperature cycling. 

2. Model, data and methods 

For the reactive system in a single rock pore we employ a spatially homogeneous flow model; 
in other words, a continuous stirred tank reactor (CSTR) paradigm. An a posteriori assessment of 
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X — Y 
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Figure 1. The quadmplex RNA species Z. The duplex replicator species are XY and xy, which for the purposes of 
the simulation are equivalent. See ref. [3 | for more details. (Online version in colour.) 



the rationale and validity of this model is given in section |4jl. The following dynamical system 
models the coupled evolution of the reactant concentrations in (Rl) and the temperature: 



V 



dc v 
~dT 



V- 



.dc w 



VC 



dt 
dT 

~dt 



-\/A;i(T)c v c w + F(cyf - c v ) 
-V2ki{T)c yCw + F(c 

w,f ) 

(-AHjVhWcyC* 

- FC(T - T t ) - L(T - T a 



(1) 
(2) 

(3) 



where the reaction rate constant fci(T) = A\ exp(—Ei/(RT)), c v is the concentration of thiosul- 
fate and c w is the concentration of H2O2. The symbols and notation used in these and following 
equations are defined in table [2j 

We used the following minimum subset of reactions from ref. [3] that can produce duplex 
RNA by autocatalysis: 



X + Y + XY 



k 2 . 



Z 2XY, 



(R2) 
(R3) 



where X and Y represent RNA complementary single- strands, the RNA duplex XY is the repli- 
cator and Z is the RNA quadruplex, the form of which is cartooned in figure [TJ In the reaction 
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volume these species evolve as 



dc 

V-gt = -\/2fc 2 (T)c x c y c xy + F(c x , f+y , f - c x+y ) (4) 
dc 

V-^ = -Vk 2 (T)c x c y c xy + V2k 3 (T)c z + F(c f , xy - c xy ) (5) 



dc 

V ~dJ = Vk ^ c * c y c *y ~ vk ^ T ) c * ~ ( 6 ) 



where the notations c x+y and c Xj f +y ,f are shorthand for c x +c y and c Xj f+c yi f, k 2 (T) = A 2 exp(— E 2 /(RT)) 
and k 3 (T) = A 3 exp(—E 3 /(RT)). (These are to be considered as empirical rate constants; (R2) 
as written is not intended to imply that an elementary three -body collision occurs.) 
The enthalpy balance includes contributions from all reactions: 

VC^ = (-AH^VhiT)^ 

+ (-AH 2 )Vk 2 (T)c x c y c xy + (-AH 3 )Vk 3 (T)c z 

— FC(T — T f ) — L(T — T a ). (7) 

We derived the activation energies, pre-exponential factors, and reaction enthalpies from the rate 
data in Q, see the Appendix. Equations ([T|)-(|3]) (no RNA supply) and equations ([l])-([2]) and 
Q-Q (RNA single- strands and duplexes supplied), were integrated using a stiff integrator from 
reasonable initial conditions. The stability of steady state solutions was assessed over a range of 
the inflow temperature, Tf, by solving the corresponding eigenvalue problem and flagging points 
where an eigenvalue changed sign. These bifurcation points then were followed over a range of 
the thermal conductance L to obtain stability maps. Numerical values of the fixed parameters are 
given in tableland are discussed further in the Appendix. 



3. Results 

First we examined the behavior of the standalone H 2 02/thiosulfate system. Figure [2|a) shows 
the computed stability map. The regime of threefold multiplicity mapped by the locus of saddle- 
node bifurcations is included for completeness, but this regime is largely irrelevant for our pur- 
poses since the temperature is too low to drive RNA replication. At each point within the Hopf 
bifurcation loop the stable solution is a limit cycle. The time series for a selected point in the loop 
is shown in figure |2jb). 

When reactions (R2) and (R3) are included, the stability map and temperature time series are 
indistinguishable numerically from those in figure [2] because the contributions of the second and 
third terms on the right hand side of equation ([7]) are relatively insignificant. 

However, time series data rendered in figure[3]for the concentrations of single-strands, duplex, 
and quadruplex show that the f^CVthiosulfate thermochemical oscillator can indeed drive the 
RNA replication process. This time series was computed for a slightly lower value of L and has a 
shorter period and higher temperature amplitude than that in figure |2jb). 

The period is determined by the ratio 4> = (CE 1 )/(Rc Wtf (—AH{)), and in general a large 
value of the numerator (high specific heat and activation energy) corresponds to longer cycling 
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Figure 2. Computed data from equations Q-Q: (a) Hopf (solid line) and saddle-node (dashed line) bifurcation 
loci; (b) Time series for Tf = 282 K, L — 2.4 mW; the oscillation period is 108.7 s. (Online version in colour.) 

times and a large value of the denominator (high specific reaction enthalpy) corresponds to shorter 
period oscillations. The physical parameters V/F (mean residence time), Tf, T a , and L also affect 
the period, as might be inferred from the locus of Hopf bifurcations in figure |2]^a) 112611 . The 
hydrogen peroxide oscillator has a variable period of the right order — around 80 to 110 seconds 
— to drive the replication of small RNAs. If the period is too long the RNAs may decay faster 
than replication can amplify them. If the period is too short the strands do not separate completely 
and replication may fail. 

But the period is not tuned to the putative doubling time for replicant concentration as was 
the thermal cycling time in [3]. As our system is an open flow system operating at dynamic 
steady state, where it has 'forgotten' the initial conditions, the doubling time is a different quantity 
from that for the closed batch system in [3]. Here it is defined as a time interval during which 
the cumulative replicant concentration at the outlet, c xy + c z , becomes equal to twice the inflow 
concentration c xy t . In equations ([T]), ([2]) and <(4|) — <(V]) the variables are evolved self-consistently and 
the doubling time is temperature-dependent, and therefore changes continuously. We can make 
use of average quantities though. 

From computed datafile, the replicant concentration c xy +c z integrated over one period of 99.3 s 
(using the trapezoidal rule on 50,000 data points for high accuracy) is 2.281 x 1CT 5 mol s/1, which 
yields an instantaneous average concentration c xy + c z = 2.297 x 10~ 7 mol/l and an average 
doubling time of 0.31 s, for c xy f = 3.5 x 1CT 8 M. This does not infer that replication is exponential, 
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Figure 3. Time series from equations ([T}, Q and for Tf = 282 K, L = 2.3 mW; period of the oscillation is 

99.3 s. (Online version in colour.) 



it simply means that on average it takes 0.31 s for the cumulative amount of replicant produced to 
become equal to twice the feed amount. 

Another measure of the system's performance is the ratio c xy + c z /c xy f = 6.56. This means 
the system has been designed so that in steady flow state, for the ideal case where no decay of 
replicant occurs within the reaction volume, the mean outlet replicant amount is 6.56 times the 
inlet duplex amount. 

Let us consider the phase relationships in figure [3] over one period. Single- strands are con- 
sumed as the temperature rises because duplex production by (R3), which requires heat, increases. 
Correspondingly the concentration of quadruplex falls, almost to zero near the temperature max- 
imum. Single-strands begin to accumulate, but there is only a small bump before quadruplex 
production picks up as the temperature declines. The presence of this bump means that the os- 
cillation is actually quasiperiodic, and reminds us that the dynamical system of equations ([l])-(|2]) 
and (]4])-(|7]> is capable of complex periodic and even chaotic behaviour. This may lend additional, 
powerful capabilities to a molecular replicating system. For example, biperiodic temperature re- 
sponse is capable of replicating two different RNA species, and nature may well have done exactly 
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Figure 4. Onset of oscillations via slow drift of parameters, (a) p 
and p2 = T. d . (Online version in colour.) 
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that in the primordial rock pores. Duplex concentration attains a broad maximum, then declines in 
phase with the temperature as reaction (R2) is favoured. The inflow delivers more single-strands, 
which accumulate even though they allow more production of quadruplex. 

3.1. Development of self-sustained oscillations 

We have addressed also the question of how self-sustained H 2 0 2 thermochemical oscillations 
may have arisen spontaneously in rock pores on the early earth. The fuzzy answer is that param- 
eters may drift quasistatically into an oscillatory regime, such as the loop of Hopf bifurcations 
in figure [2ja). A precise answer is provided by simulating various scenarios where one or more 
tunable parameters p in equations [l]-[3] is given the following increasing time dependence, such 
that p approaches a constant value: 

p = poo (1 - a exp (-7*)) . (8) 

A time series from equations (jT]) — ((3]) using ([8]) for one or more parameters is expected to show 
the temperature increasing slowly, then beginning to wobble, then settling into self-sustained os- 
cillations. Some examples are shown in figure |4} Of course, such evolutions could proceed for 
millions of years before the onset of periodic behaviour, but we we have set the time constant 7 in 
equation ([8]) more conveniently. 



4. Discussion 

4.1. Validity and relevance ofCSTR paradigm 

The CSTR paradigm described by equations (jT]) — ((2]) and (|4])-(|7]) can be thought of in two ways. 

On the one hand it provides a recipe for a laboratory experiment in a microreactor cell. The 
flow rate F, inflow concentrations c v f, c Wj f, c x f +y f and c xy f, inflow temperature T f and wall or 
ambient temperature T a are tunable by the experimenter, the wall thermal conductance L and 
volume V are set by design, and thermokinetic parameters and reaction enthalpies of the RNA 
reactions can be manipulated by engineering the ribonucleotide sequences and strand lengths. A 
variety of continuous-flow microreactors are already in use for various calorimetric applications 
in biotechnology, and the field is developing rapidly ll27l . 

On the other hand we need to address the adequacy of the CSTR as a model for reactions in 
a porous rock. In an actual physical situation we assume a flow through a porous rock structure 
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comprising a large array of individual pores, having patterns of connection determined by the 
character of the rock structure. Provided that the system is sufficiently ergodic, and input flow is 
steady in volumetric velocity and chemical composition, it is reasonable to assume that the time 
behaviour in a single pore as reported above gives results at least qualitatively representative of 
the larger array over times longer than the period of the thermochemical oscillator. 

Further, the use of a CSTR model for a pore implicitly assumes a typical residence time 3> 
the mixing time. Mixing is achieved mainly by convection, and in this respect thermoconvection 
times of 3-7.5 s were reported by Mast et al. Il2~8~1l . The mean residence time in our system is 20 s. 
The perfect mixing assumption is considered to hold adequately when the residence time is around 
5-10 x the mixing time, so the CSTR model can be a reasonably good approximation. The 15 s 
thermoconvection time reported by Braun and Libchaber [|29l is too long for the perfect mixing 
assumption to hold, and in that case the appropriate model is a system of reaction-convection- 
diffusion equations and boundary conditions. In that case reactive and convective thermal oscil- 
lations would couple in important and interesting ways, but that is for future study. We prefer, at 
this stage, to avoid the sheer 'tyranny of numbers' that still plagues the numerical analysis of such 
systems, and isolate the reactive thermal oscillations by studying a system for which the perfect 
mixing assumption holds. 

In experiments the mixing behaviour of the system would be characterised by using a tracer 
to measure the residence time distribution, or or by determining the period of time necessary for 
the system to achieve a desired level of homogeneity. We note that Imai et al. 11301 confirmed 
experimentally the well-stirred condition for a flow microreactor that simulated a hydrothermal 
environment, in which elongation of oligopeptides was achieved at high temperature and pressure. 

4.2. RNA stability to hydrogen peroxide 

The response of RNA to H2O2 is a mixed story. Some small RNAs have been found to be 
very stable; for example a 109-nucleotide RNA which is induced by oxidative stress in E. coli has 
an in vivo half-life of 12-30 minutes [|3T| . There is evidence for modern RNA dysfunction, but 
not degradation, caused by oxidation initiated by H 2 0 2 Il3~2l . In [|33l it was found that small RNA 
(~50 bases) of yeast is stable to H2O2, but pyrite is reactive in RNA degradation and hydroxyl 
(OH) radical generation. 

Much more is known about the effects of H 2 0 2 on DNA, in the context of disease caused by 
free radical damage. In [34] it was found that 2'deoxyguanosine was hydroxylated at C-8 when 
included in the thiosulfate-hydrogen peroxide system. But it is an interesting fact that H 2 0 2 does 
not react directly with DNA, and DNA is not damaged in its presence unless transition metal ions 
are present 11351 . In the presence of transition metal ions OH radicals are generated by the Fenton 
reaction: Fe 2+ + H2O2 — > Fe 3+ + OH + OH . However, the high and indiscriminate reactivity of 
the OH radical limits its ability to damage biomolecules, because it is more likely to be scavenged 
by the iron before attacking DNA or RNA: OH + Fe 2+ ->■ Fe 3+ + OH . 

Oxidized RNA evidently may still replicate. For example, oxidized mRNA has been success- 
fully converted to cDNA by reverse transcriptase, with mutations induced in the cDNA by the 
oxidized RNA bases 136). 

In summary, although H 2 0 2 may damage or modify RNA (as could many other possible in- 
gredients of primordial soup), there is no reason to suppose that it cannot replicate and pass out 
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of the reaction zone faster than it is damaged or degraded. RNA that is modified by the action of 
H 2 0 2 in such a way that confers resilience to H 2 0 2 damage would, of course, be selected for. 

4.3. Occurrence and supply of substrate RNA oligomers 

Replication of RNA duplexes by complementary strand pairing raises the thorny question of 
how such strands could have been produced in the prebiotic primordial soup. It has been shown 
that 200-mers of RNA can be polymerized and concentrated in a thermal gradient [37|, but this 
does not produce complementary strands preferentially. Polynucleotides > 50-mers can be syn- 
thesized on montmorillonite [f38l . which also facilitates homochiral selection ||39l , and this would 
seem a stronger hypothesis. 

A supply of complementary template oligomers could be provided by a variation on the surface- 
promoted replication procedure in [|40l , when competitive pathways are included for consumption 
and production of complementary free RNA species X and Y as indicated in figure [5J Here, 
too, replication may be driven by an H 2 0 2 thermochemical oscillator, because strand separation 
requires heating and annealing, ligation, and immobilization of single strands is promoted by cool- 
ing. 



Figure 5. Scheme for supply of template oligomers. The surface-bound single polynucleotide strands S-X and S-Y 
bind n complementary nucleotide fragments N and N' from solution, which are ligated irreversibly to form S-XY 
and S-YX. Heating promotes release of free X and Y and slower release of duplex XY from S, and melting of free 
XY to X and Y. The free species can participate in the duplex replication scheme (R2) and (R3). 

In another pre-biotic scenario, lipid molecules in a dehydrated environment promoted conden- 
sation of nucleic acid monomers, and an alternating wet phase replenished the monomers PTI . 
Such a system also would be enhanced by the H 2 0 2 thermochemical oscillator, which would pro- 
vide periodic heating to drive the dehydrations and cooling to allow rehydrations. 

4.4. Effects of other reactions involving H2O2 

Hydrogen peroxide may undergo other competing reactions; for example, the Fenton reaction 
in the presence of Fe 2+ /Fe 3+ discussed above, which is also exothermic. In this case there would 
be an extra reaction rate term in the species balance for hydrogen peroxide and in the thermal 
balance. Such a system has an enhanced propensity towards quasiperiodic oscillations. 

It is noteworthy that the Cu(II)-catalysed H 2 0 2 /thiosulfate ion reaction displays pH oscillations 
isothermally Il42ll . implying that the chemistry alone is sufficiently nonlinear to permit oscillatory 
dynamics, and in any case more complex than (Rl) suggests. Spontaneous pH cycling could well 
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be an assistant driver for RNA replication as it would provide scope for favourable conformation 
changes. 

In summary, we have proposed and explored a previously unrecognized mechanism for driving 
a replicating system on the prebiotic earth. The chances of a hydrogen peroxide thermochemical 
oscillator arising spontaneously on the early earth in the presence of nucleotide precursors are 
perhaps very small. However, the earth is large, rock pores are innumerable, and there was plenty 
of time on the early earth for improbable events to happen. 

Braun and Libchaber ll29ll are strong proponents of an interdisciplinary approach between bio- 
chemistry and geophysics for understanding the origin of life at the molecular level. We have 
added some insights that originated in chemical engineering, showing that an interdisciplinary 
approach is indeed fruitful. 



Table 2. Nomenclature and numerical values. 



Symbol (units), definition 


Quantity 




Values 






no RNA 


with RNA 


A (M n_1 s _1 ), pre-exponential factor 


At 


1.63 x 10 


LU 1.63 x 10 1U 




A 2 


n/a 


1.69 x 10 12 




A 3 


n/a 


1.0 x 10 39 


c (M), concentration 


C v ,f 


0.9 


0.9 




Cw,f 


1.35 


1.35 




Cx,f+y,f 


0 


5.0 x 10 -7 




Cxy.f 


0 


3.5 x 10 -8 


C (JK -1 1 -1 ), volumetric specific heat 


c 


3400 


3400 


E (kJ mol -1 ), activation energies 




68.12 


68.12 




E-2 


n/a 


-10 




E, 


n/a 


260.14 


F (/ils -1 ), flow rate 


F 


0.50 


0.50 


AH (kJ mol -1 ), reaction enthalpies 


AH! 


-612.5 


-612.5 




AH 2 


n/a 


-169.6 




AH 3 


n/a 


260.14 


L (WK _1 ), wall thermal conductance 








n, sum of species reaction orders 








T (K), temperature 


% 


283 


283 


V (fiY) pore volume 


V 


10.0 


10.0 



Subscripts 



1,2,3 pertaining to (R1),(R2),(R3) 

a of the ambient or wall temperature 

f of the flow 

v, w of thiosulfate, H2O2 

x, y, z of single-strands, duplex, quadruplex 
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Appendix: Choices of numerical values 

• The weighted average volumetric specific heat C = 3400 J KrH -1 is considerably lower 
than that of pure water. The volumetric specific heat of sea water is about 3850 J K _1 l _1 and 
that of H 2 0 2 is about 2620 J K -1 1 _1 , from which we arrive at the given value of 3400 J K _1 l _1 

• The activation energies and pre-exponential factors of the RNA reactions were obtained 
from the data in 01 using 

l/7\ - 1/T 2 

then A = kexp (E/RT). The activation energy for the association reaction (R2) is zero 
in (3]| but we used the small negative value of lOkJ/mol in the calculations to reflect the 
fact that the rates of association reactions of biological macromolecules often decrease with 
temperature, because due to thermal motion a smaller faction of energetically favourable 
collisions result in reaction. The data in [|3l gave a pre-exponential factor for reaction (R3) 
in the main article of 10 45 s _1 . On the basis of only two data points this should not be 
taken too literally; in practice we had to reduce it to 10 39 s _1 to couple the reaction to the 
thermochemical oscillator. 



Acknowledgement: Rowena Ball is recipient of Australian Research Council Future Fellowship 
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Figure captions and short title for page headings 

Figure 1. The quadruplex RNA species Z. The duplex replicator species are XY and xy, which 
for the purposes of the simulation are equivalent. See ref. 0 for more details. (Online 
vesion in colour.) 

Figure 2. Computed data from equations ([T])-([3]): (a) Hopf (solid line) and saddle-node (dashed 
line) bifurcation loci; (b) Time series for T f = 282 K, L = 2.4 mW; the oscillation period is 
108.7 s. (Online version in colour.) 

Figure 3. Time series from equations {j}, ^ and (|4])-0 for T f = 282 K, L = 2.3 mW; period of 
the oscillation is 99.3 s. (Online version in colour.) 

Figure 4. Onset of oscillations via slow drift of parameters, (a) p = c Wi f, (b) p = T a , (c) p = F, 
(d) both pi = c W(f and p 2 = T a . (Online version in colour.) 

Figure 5. Scheme for supply of template oligomers. The surface-bound single polynucleotide 
strands S-X and S-Y bind n complementary nucleotide fragments N and N' from solution, 
which are ligated irreversibly to form S-XY and S-YX. Heating promotes release of free 
X and Y and slower release of duplex XY from S, and melting of free XY to X and Y. 
The free species can participate in the duplex replication scheme (R2) and (R3). 

Short title: Primordial RNA replication 
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